arXiv: 1509.07112v2 [quant-ph] 7 Dec 2015 


Quantum Inf Process manuscript No. 

(will be inserted by the editor) 


Quantum Walk on the Line through Potential Barriers 

Thomas G. Wong 


Received: date / Accepted: date 


Abstract Quantum walks are well-known for their ballistic dispersion, travel¬ 
ing 0(t) away in t steps, which is quadratically faster than a classical random 
walk’s diffusive spreading. In physical implementations of the walk, however, 
the particle may need to tunnel through a potential barrier to hop, and a naive 
calculation suggests this could eliminate the ballistic transport. We show by 
explicit calculation, however, that such a loss does not occur. Rather, the 0(<) 
dispersion is retained, with only the coefficient changing, which additionally 
gives a way to detect and quantify the hopping errors in experiments. 
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1 Introduction 

Quantum walks are the quantum analogues of classical random walks. In the 
same way that classical random walks, or Markov chains, serve as useful al¬ 
gorithmic tools for classical computing [T], quantum walks are the basis for 
several important quantum algorithms, including search 03 ], element dis¬ 
tinctness [1] , and boolean formula evaluation For a comprehensive review 
of quantum walks, see [7]. 

Before quantum walks were applied algorithmically, however, their basic 
properties were investigated, beginning with how they evolve on the one¬ 
dimensional infinite lattice. This began with Meyer [5] in the context of quan¬ 
tum cellular automata, who showed that a discrete-time quantum walk neces¬ 
sitates an internal degree of freedom to evolve non-trivially. Interpreting this 
internal degree of freedom as spin, he showed that the evolution’s continuum 
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limit is the Dirac equation of a spin-1/2 particle. Later, this work was trans¬ 
lated into the language of quantum walks and extended by Ambainis et al. [9l 
[in], revealing an important distinction between classical and quantum walks: 
quantum walks spread quadratically faster than classical ones. This was also 
shown to be true for continuous-time quantum walks m , which do not require 
the internal degree of freedom [T2|. Since these investigations, there has been 
continued interest in one-dimensional quantum walks; just some examples in¬ 
clude [HdllllllllllTlIIHI. 

Spurred by this substantial theoretical interest, several experimental groups 
have implemented one-dimensional quantum walks in a variety of physical 
systems, including linear optical resonators m, nuclear magnetic resonance 
samples [20] , photons in waveguide lattices iniEa, photons in interferometric 
networks [53], optically trapped atoms [53], and linearly trapped ions HSIIMI 
HZ]. For most of these, the walk occurred in discrete-time. 

In this paper, we also focus on the discrete-time quantum walk on the one¬ 
dimensional infinite lattice, but with two important differences that we will 
introduce below. To begin, we use computational basis states labeled by the 
integers {..., |—2), | —1), |0), |1), |2),... } for the position of the particle. Since 
there are two possible directions for the particle to hop (left and right), we 
take the particle to be a two-component spinor with computational basis states 
11^), |—>■)}. Then the state of the system is a vector in C°° 0 C^, where the 
first tensor factor corresponds to the position of the particle, and the second 
to its internal spin or “coin” state. The quantum walk is given by repeated 
applications of the unitary operator 

5(7 0 if), 

where H is the “Hadamard coin,” which operates on the coin space and takes 

H\'^) = ( 1 ^) + 1 ^)) 

H\^) = ( 1 ^) - |->)), 

and S is the “flip-flop” shift 

= I* ~ 1, ~^) 

‘5'lb ~^) = I* + 1, ^)> 

which flips the internal state of the walker when hopping. The use of the flip- 
flop shift S is the first important distinction of our analysis, which contrasts 
with typical investigations of ID quantum walks that instead use the “moving” 
shift, where the direction of the walker is unchanged after the hop. There are 
several reasons for our choice of the flip-flop shift. First, it is the natural 
operator for a quantum walk using a macroscopic Bose-Einstein condensate 
[28| . where the condensate atoms shift by exchanging photons with two laser 
fields, undergoing a stimulated Raman transition. This imparts a momentum 
kick to the atoms that physically shifts them while flipping their internal states. 
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Lattice Position 


Fig. 1 Probability at even lattice sites for an initially left-moving particle localized at 0 
after 100 steps of the quantum walk. The solid black curve is with the flip-flop shift, and 
the dashed red curve is with the moving shift. Note the probability at odd sites is zero. 


A second reason for using the flip-flop shift is that it has not been previously 
analyzed for the ID quantum walk. A third reason is that the flip-flop shift is 
Hermitian, which will be important in a moment when we introduce potential 
barriers. A fourth reason for using the flip-flop shift is algorithmic: it allows 
quantum walks to search graphs quickly, whereas the moving shift, by contrast, 
yields a slow search algorithm [53]. Furthermore, as a fifth reason, defining 
the flip-flop shift on non-Euclidian lattices, such as the complete graph |3D| or 
hypercube |5], is natural, whereas it is unclear what the moving shift would be. 
Finally, a sixth reason is that for a quantum walk alone (i.e., simply applying 
the coin and shift operator, so no oracle or algorithmic modifications), the flip- 
flop shift yields the same probability distribution as the moving shift, apart 
from a directional flip. 

To illustrate this final point, say the quantum walk starts in the middle of 
the lattice (he., at position 0), pointing left: 

IV’o) = |0) 0 1^). (1) 


The probability distribution of the quantum walk after 100 steps, for both 
the flip-flop and moving shifts, is shown in Figure]^ Surprising, the flip-flop 
shift causes the walk to become right-moving, while the moving shift retains 
the initial left-moving behavior. In addition, the probability distributions are 
exact reflection of each other. To see why this occurs, note that after two steps 
of the quantum walk, the state of the system with each of the shifts is 


flip-flop) 
|'02,moving) 


|-2)®K)-|0)®(M-K))-|2)®|^) 

|—2) (g) 1^) -I- |0) (g (!■<—) -I- |—>■)) — |2) (g |—>•) 


At this point, the probability distributions of the walks are identical. It is 
during the third step that the directional flip occurs. To see this, focus on 
the coin state at the middle vertex |0), which for the flip-flop shift is (!•(—) — 
|—>'))/'\/2, and for the moving shift is (!•(—) + \^))/V2. When the Hadamard 
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Lattice Position 


Fig. 2 Probability at even lattice sites after 100 steps of the quantum walk (solid black) 
and classical random walk (dashed red) with an unbiased initial state in the middle of 
the lattice (so the quantum starts in and the classical starts with probability 1 at the 
middle). Note the probability at odd sites is zero. 


transform acts on the flip-flop shift’s middle coin state, it becomes |—>■), so 
the walk has become right-moving. For the rest of the walk, it continues to be 
right-moving, with the probability distribution mirroring the moving shift’s 
left-moving walk. When it acts on the flip-flop shift’s coin state, however, it 
becomes |^), so it retains its left-moving bias. 

Since the quantum walk is a purely real operator, if the system starts the 
unbiased initial state localized in the middle of the lattice, 

|V>o) = |0)®^(|^)+*K)), (2) 

then the flip-flop and moving shifts will yield exactly the same probability 
distribution. This distribution after 100 steps of the walk is shown in Figure 
along with the classical random walk’s probability distribution; this reveals 
the quantum walk’s linear, ballistic propagation, which is quadratically faster 
than the classical walk’s square-root diffusion. That is, after t steps, the quan¬ 
tum walk spreads an average distance 0(t) while the classical random walk 
spreads 0(-\/t). Note since the quantum walk will not mix real and imaginary 
components, rather than using the unbiased initial state (§ , it suffices to an¬ 
alyze the left-moving initial state Q instead, and we assume this for the rest 
of the paper. 

Rather than just analyzing this ID quantum walk, we augment it by in¬ 
troducing potential barriers that the quantum walker must tunnel through in 
order to shift mm- That is, we replace the flip-flop shift with 

S^aS + /3/, 

so there is some amplitude a that the particle successfully tunnels, and some 
(related) amplitude /3 that it fails and stays put. As previously noted, the flip- 
flop shift S is Hermitian, so the modified shift is unitary provided |ap-|-|/3p = 1 
and a/3* -I- /3a* = 0. If the moving shift was used instead, then this operator 
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Fig. 3 One step of the quantum walk with potential barriers (Le., the Hadamard coin 
followed by a shift/identity) for a (a) left-moving and (b) right-moving particle. 



Fig. 4 Probability at lattice sites for an initially left-moving particle localized at 0 after 
100 steps of the quantum walk. The solid black curve is the barrier-free (q = 1 and ^ = 0) 
walk (with only even positions shown), and the dashed red curve is with a = cos(0.8) and 
(3 = 2sin{0.8) (with all positions shown). 


would not be unitary, and would hence not be an quantum operator. With the 
potential barriers, a complete step of the quantum walk is 

U ={aS + pi){I®H), (3) 

and its action on a left- or right-moving particle is shown in Figure]^ Modi¬ 
fying the shift is a second distinction of this paper from previous work on ID 
quantum walks. 

There are other physical processes where this modihed shift can arise. For 
example, a quantum walk using atoms trapped in an optical lattice |24| can 
be shifted by translating the trapping potential. Even if the trapping potential 
is moved at constant speed, its acceleration is nonzero at the beginning and 
end of the motion [^, and this acceleration effectively causes the trapping 
potential to tilt [34], resulting in the particles staying put. As another example, 
a quantum walk using a macroscopic Bose-Einstein condensate |28| shifts when 
a stimulated Raman transition imparts a momentum kick to the atoms. The 
uncertainty as to whether all the atoms have shifted once can cause a portion 
of the condensate to stay put. 

In Figure]^ we plot the probability distribution of the quantum walk after 
100 steps both with and without the potential barriers. Recall we are assuming 
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for the rest of the paper that the particle is initially left-moving and localized 
in the middle of the lattice according to so the particle becomes right- 
moving. This figure shows that the potential barriers slow down the walk, but 
how much? As a naive calculation, consider the parameterization a = cos (p 
and /3 = isint/j with p small, so a « 1 and /3 « ip. Then from Section 4.5.3 
of Nielsen and Chuang [35], the error between the perfect shift S and the 
one with potential barriers (aS + j3I) is roughly p, and after t steps of the 
quantum walk operator @ , the error is upper bounded by pt. Since this bound 
allows the error to increase as the system evolves, one might expect that the 
characteristic behavior of the quantum walk deteriorates, perhaps also loosing 
the ballistic propagation. 

In the next section, however, we explicitly show that this does not happen. 
We analytically determine the evolution of the particle and show that the 0(<) 
dispersion is retained so long as the potential barriers do not asymptotically 
cause the quantum walk to stop hopping. That is, the potential barriers do not 
eliminate the characteristic ballistic transport of quantum walks, only scaling 
it to some degree. Furthermore, this scaling gives a way to detect and quantify 
if these errors are occurring in experiments. Finally, we end with a comment 
about continuous-time quantum walks similarly retaining the 0(t) dispersion 
with potential barriers. 


2 Analysis with Potential Barriers 


Our analysis follows the Fourier method of wm, and we begin by writing the 
wave function at position n and time t as a two-component spinor 




ypR{n,t)J ’ 


where L and R denote the left- and right-moving coin components, respec¬ 
tively. Using Figure]^ one step of the quantum walk operator ^ causes this 
state to become: 

V'(n,t-bl) = ^ il)(n+l,t)+ 

Let us call these matrices M+, M_, and Mq, where the subscript respectively 
denotes if the position shifts right, left, or stays put. Then 


ip{n, t -f 1) = M+tf}{n — 1, t) -I- M-ip{n -I-1, t) -I- Motp{n, t). 


Now consider the Fourier transform of the position wavefunction, which gives 
the momentum wavefunction tp with k € [—7r,7r]: 
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Then the momentum wavefunction at the next timestep is 

t + 1) = ^ ip{n, t + 

n 

= {M^'ip(n — 1, t) + + 1, t) + Mo'ip(n, t)) e 

n 

= (M+e*'= + + Mq) ^ t)e 

n 

= Mk^{k,t), 

where we have defined 

Mk — Tf_|_e M—6 Mq — — ( —ik o —ik o 1 ■ 

+ \/2 + /3 ae - (3) 

Using this recursive relation, the momentum wavefunction at later time is 
related to the momentum wavefunction at initial time by 

= {MkY'ijjik,0), 

where the initially left-moving particle at the center of the lattice Q has 
momentum wavefunction 

^(^>0) = (^Q 

for all k. So to find the evolution, we diagonalize Mk, which makes it easy to 
take powers of it. To find its eigenvalues A, we solve the characteristic equation 

0 = det(Mfe — A/) = A^ — aV2cos{k)\ 1. 

This yields eigenvalues 

, a cos k ± i\/2 — cP- cos^ k ., 

A . --U 


where 


V2 


a cos k . V2 — A cos^ k 

cosojk = ^ , smwfc = - 


(4) 


Now let us find the eigenvectors of Mk, which satisfy 


Mk 




The top row of this is 


^ [(ae*'" + li)u+ (-ae*'" + 0)v\ = e 

and with some algebra, it becomes 
a-k/3e-*'= - 




a — Pe 


— ik 


-U = V. 
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Thus the (unnormalized) eigenvectors of are 




Let us normalize this by finding the square of its norm: 


= UU +VV 


= 1 + 


+ \ f a* + j3*e^^ - 


a — j3e 


a* — P*e'^ 


Now we assume that a is real and /3 is pure imaginary; for example, they 
could be parameterized as a = cos(/3 and f3 = isint^. Since lap + |/3p = 1, 
this implies that a^ — /3^ = 1. So the square of the norm of the eigenvectors 
becomes 


^ 2 + 1 — 2iaj3 sin(A:) — 2-\/2a cos (±Wfc — k) + 2-\/2z/3 sin 

1 + 2za/3 sin(fc) 


(±Wfc) 


4 — 2\/2a cos (±Wfc — fc) ± 2\/2ij5 sin (wfe) 
1 + 2za/3 sin(A:) 

Thus the normalized eigenvectors of Mk are 


l^±) = 


I 1 + 2za/3 sin(fc) 

4- 2v^a cos (±Wfc — fc) ± 2\/2ij3 sin (wfc) 


Q:— 


Using these eigenvalues and eigenvectors, we diagonalize Mk and take matrix 
powers of it: 

(Mfc)‘ = (A+)‘|V'+)(V'+| + (A_)‘|z/:_)(z^_|. 

Then acting by this on the initial momentum wavefunction 0) = (1,0)''', 
we get the momentum wavefunction at time t: 


t) = e*' 


, 1 + 2za/3 sin(A:) 

4 — 2'j2a cos {tOk — k) -\- 2-\/2z/3 sin (wfc) 


Q:— 


, 1 + 2za/3 sin(fc) 

4 — 2'j2a cos (wfc + fc) — 2\/2i/3 sin (ojk) 


From this, we identify the left and right components: 


— ifc 


^L{k,t) = e*' 


( 1 + 2za/? sin(fc) 

4 — 2y/2a cos {ujk — k) + 2\/2ij3 sin (wfc) 

_ l + 2za/?sin(fc) _ 

4 — 2\f2ot cos (wfc + k) — 2-\/2z/3 sin (wfc) 
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i’Rik.t) 


1 + 2ia(3 sm{k) 




4 — 2^/2a cos (w^ — k) + 2\/2ij3 sin (wfc) o — /3e“*^ 

1 + 2za/3 sin(fc) a + /3e“*^ — 

4 — 2\/2q; cos {uk + k) — 2y/2il3 sin (wfe) a — /3e“*^ 

These are the left and right components of the momentum space wavefunc- 
tions. 

To go to position space, we take the inverse Fourier transform of the mo¬ 
mentum wavefunction: 


ip{n,t) = 


dk 

27r 


^{k,t) 


—ikn 


Let us begin with the left direction. The inverse Fourier transform is 




1 -|- 2ia/3 sin(fc) 


^—ikn 


27r 4 — 2-\/2acos (wfc — fc)-I-2-\/2i/3sin (wfe) 

r 12m/3 sin(fc) 

2tt 4 - 2V2a cos (wfe + k) — 2^/2il3 sin (wfe) 

Let’s focus on the first integral: 


„—ikn 


dk 


,iuikt 


1 -|- 2iaj3sm{k) 


—e 


— ikn 


J-TT 27r 4 — 2\/2a cos — A:) -I- 2\/2i^ sin (w^) 

Let A; = u -I- TT. Then this integral becomes 


du 


1 -I- 2*a/3 sin(M -|- tt) 


27r 4 - 2-sf2a cos - [u + tt)) -|- 2v^i/3 sin (uJu+tv) 

where from 

cos(tt -I- tt) — cos u 

COS ■\/^ 

implies that 

a;„+^ = cos“^(-cosa;„) =7r-a;„. 

So this integral simplifies to 


^ — i{u+Tr)n 


r _ 1 - 2ial3 sin{u) _^ 

/- 27 I- 27r 4 — 2\/2a cos (a;„ -I- u) -b 2'/2ij5 sin (w„) 




Since the integrand is periodic with period 27r, we can shift the bounds to 


1)"+* 


1 — 2ia/3sin(ii) 


27r^ 4-2\/2 a cos (wu -b u) -b 2y/2iP sin (a;„) 
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Lattice Position 


Fig. 5 Probability at each lattice site for an initially left-moving particle localized at 0 after 
100 steps of the quantum walk with a = cos(0.8) and ^ = 2sin(0.8). The solid black curve is 
from iteratively simulating the evolution, and the dashed red curve is from integrating © 
and §. The two methods are in agreement, so the curves overlap. 


Combining this with the second integral, the position-space wavefunction in 
the left direction is 




+ 


!)»+« l-2m/3sin(fc) ^-i(u^a+kn) 

4 - 2V2a cos (wfe -I- /c) -I- 2\/2if3 sin (wfc) 
r _ l + 2ial3s,m{k) _ 

J_^ 27r 4 - 2y/2a cos (wfc + k) — 2^/2ij3 sin (wfc) 


Similarly, the right-pointing wavefunction is 


fpR{n,t) 


_1 - 2ial3 sin(fc)_ 

4:-2V2a cos {ujk + k) + 2\/2if3 sin {ujk) 


dk 


a + 

1 -I- 2iaf3 sin(fc) 


l_^ 27r 4 — 2^/2a cos {uk + k) — 2\/2il3 sin (wfe) 


a — j3e 


( 6 ) 


Although these formulas ([^ and are complicated, they can be numerically 
integrated and verified to agree with a numerical simulation of the quantum 
walk; this is shown in Figure where the two probability distributions gen¬ 
erated by each method overlap in agreement. Also note that the wavefunction 
at every other lattice site is no longer equal to zero with the potential barriers, 
even though they were without the barriers. 

Now that we have the position-space wavefunctions in each direction, we 
can approximate them for large times (t 3> 1). As described in |10j . this can 
be done using the method of steepest descent, where we examine the time- 
oscillatory portion of (§ and § , which is 

^-i(ujkt+kn) _ ^-i<p(k,a)t 
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where we have defined 

, / 7 \ 

(t){k,a) =u}k + —■ 

From the analysis of nni, the part we care about is the movement of the peak 
in probability {e.g., the leading right peak in Figure]^. This is given by where 
(j) has first and second derivative (with respect to k) equal to zero. So let us 
calculate these derivatives, beginning with the first: 


d6 . n 

To find a;(,, we use implicit differentiation on Wj’s definition which yields 

, a sin k 


Wfc = 


\/2 — cos^ k 


Plugging this in, the first derivative of (f) is 


d(j) a sin k n 

dk \/2 — cos^ k t 


Now let us differentiate this to find the second derivative of (j): 


d‘^(j) a cos fc (2 — 0 !^) 
dk"^ (2 — cos^ 

Using these expressions for the first and second derivative of (p, we find they 
are both zero when 

TT Stt 7 11 Of 

k = and — = . 

2’ 2 t 

It is easily verified that this is a 2nd order stationary point, so d^(j)/dk^ ^ 0 
at these values. Thus the peak occurs at position 

a 

n = 

72 

at time t, which is the quantum walk’s characteristic linear dispersion 0(t). 
That is, the spreading remains ballistic, with the potential barriers only affect¬ 
ing the constant factor. Note when a = l^we get the expected potential-free 
peak at t(y/2. This is confirmed in Figures where we plot the location of the 
peak as the quantum walk operator (|^ is applied; the potential-free curve has 
a linear regression with slope nearly 1/72 « 0.707, and with a = cos(0.8), the 
slope is nearly cos(0.8)/72 « 0.493. 

This result can also be used to detect and quantify the hopping errors in 
experiments. If the movement of the quantum walk has a speed slower than 
1/72 lattice positions per unit time, then one might suspect that hopping 
errors are occurring of the form {aS + 131). Since the speed equals Qf/72, this 
in turn yields a and /? = f7l — 
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F'l^. 6 L3jttic0 position of tli 0 mSiXimiim probability ns stops of tli 0 cjimntum Avnlk aro 
takon with an initially loft-moving particlo localizod at 0. Tho solid black curvo is with a = 1 
and /3 = 0, and tho dashod rod curvo is with a. = cos(0.8) and (3 = isin{0.8). 


3 Continuous-Time Quantum Walks 

We end with a short comment on continuous-time quantum walks [12] , which 
evolve on fc-regular graphs by Schrodinger’s equation with Hamiltonian 

H = -7H, 

where 7 is the jumping rate (be., amplitude per time) of the walk, and A is the 
adjacency matrix of the graph {Aij = 1 if vertices i and j are connected and 0 
otherwise). The potential barriers would hinder the amplitude from jumping, 
say by e, and cause it to stay put instead. This effectively causes the adjacency 
matrix to be |31| 

A' = kel + (1 - e)A. 

But the multiple of the identity matrix can be dropped since it only con¬ 
tributes an overall, unobservable phase. Thus the system effectively evolves 
with Hamiltonian 

H' = -7(1 - e)A = (1 - e)H. 

But this is just a rescaling of the potential-free Hamiltonian, so the system 
evolves the same way, but with rescaled time t' = {1 — e)t. So assuming that 
the potential barrier yields a reduction in amplitude e that does not approach 
1 , the ballistic transport of the continuous-time quantum walk m is retained, 
similar to the discrete-time quantum walk we showed above. 


4 Conclusion 

We have analyzed the discrete-time quantum walk on the infinite ID lattice, 
replacing the typical moving shift with the fiip-flop shift, which is the one 
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actually used in quantum walk search algorithms. With this, we considered 
the effect of potential barriers hindering the hop, showing that as long as the 
barriers do not asymptotically cause the walk to stop, they do not change the 
0(t) ballistic transport of the quantum walk. This does change the coefficient 
of the dispersion, however, and this provides a way to detect and quantify the 
hopping errors. These results also apply to continuous-time quantum walks. 

Since our model assumes that the potential barriers are equal everywhere 
on the lattice and do not fluctuate, an open question is the effect of nonhomo- 
geneous or random barriers. This contrasts with previous work on nonhomoge- 
neous or random coins and comparing or combining the investigations 

would be particularly interesting. 

Acknowledgements Thanks to Andris Ambainis and Alexander Rivosh for useful dis¬ 
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